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Abstract 

Structural pattern recognition describes and classifies data based on the re- 
lationships of features and parts. Topological invariants, like the Euler num- 
ber, characterize the structure of objects of any dimension. Cohomology can 
provide more refined algebraic invariants to a topological space than does ho- 
mology. It assigns 'quantities' to the chains used in homology to characterize 
holes of any dimension. Graph pyramids can be used to describe subdivi- 
sions of the same object at multiple levels of detail. This paper presents 
cohomology in the context of structural pattern recognition and introduces 
an algorithm to efficiently compute representative cocycles (the basic ele- 
ments of cohomology) in 2D using a graph pyramid. An extension to obtain 
scanning and rotation invariant cocycles is given. 

Keywords: graph pyramids, representative cocycles of cohomology 
generators 
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1. Introduction 

Image analysis deals with digital images as input to pattern recognition 
systems with the purpose to extract information about their content, usually 
objects. Objects appear in images affected by transformations (e.g. rotation, 
zoom, projection) and noise. Topological features have the ability to ignore 
changes in the geometry of objects by extracting object properties invariant 
to elastic transformations. Simple topological features are for example the 
number of connected components, the number of holes, etc., while more 
refined ones, like homology and cohomology, characterize holes and their 
relations. 

An example application of topological features is topology simplifica- 
tion, an active field in geometric modeling and medical imaging where high- 
resolution surfaces are created through iso-surface extraction from volumetric 
representations, obtained by 3D photography, CT, or MRI. Iso-surfaces often 
contain many topological errors in the form of tiny handles. These nearly 
invisible artifacts hinder subsequent operations like mesh simplification, re- 
meshing, and parametrization. See, for example [ij. Another application 
is shape description and matching, where persistence and homology of a 
function defined on a shape have been successfully applied to extract shape 
features jij. 

A 2D image is the result of projecting a 3D scene into the image plane. 
Often the precise camera parameters are not known and still humans have 
no problem in correctly interpreting the displayed objects in the image. A 
3D object is surrounded by a reflecting surface which itself may split into 
several smaller but connected patches which can be characterized by their 
color, their texture or other visual properties. The visible part of the object's 
surface maps into a region of the image which shows the same adjacencies 
of patches as the original surface (because the camera sees the same side 
of the surface) although the geometry of the patches may change due to 
projection, due to camera or object movement or due to deformations of the 
object. Sometimes a collection of patches is completely surrounded by some 
other surface patches, e.g. a fancy soccer ball on a Spanish T-shirt (see Fig. 
H]). Although both the picture of the ball as well as the T-shirt may have a 
specific patch structure (stripes on shoulder and arms, logo) it is clear which 
subset of regions forms the ball and which regions belong to the remaining 
T-shirt. A simplified version could describe the pixels of the T-shirt I's and 
the ball's pixels by (see Fig. [T]). Then the ball is a hole in the T-shirt 
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a) 



b) 



Figure 1: a) Original image: ball on T-shirt of world champions 2010 in soccer; b) a 
segmentation of the original image using pyramids. 

as long as the surrounding of the ball is visible and not occluded by other 
objects. It may be highly difficult to uniquely and reliably identify any of 
the involved small patches individually under difficult geometric deformations 
while the overall arrangement of patches forming the ball is mostly invariant 
to these geometric deformations. How to segment patches into meaningful 
aggregations has been dealt with in many other segmentation methods and 
it is not the main emphasis of this paper. We therefore restrict ourselves in 
the following on binary images with the understanding that each region may 
be the collection of several subregions belonging together. 

Considering 2D binary images, an object is defined by a connected set (4- 
connectivity) of foreground pixels. A region adjacency graph (RAG) encodes 
the adjacency of regions in a partition. The holes in a RAG associated to an 
object of a 2D binary digital image can be characterized by establishing an 
equivalence between all the cycles as follows: two cycles are equivalent if one 
can be obtained from the other by joining to it one or more degenerate cycles 
(cycles with exactly 4 edges). For example, there is only one equivalence 
class for the foreground (gray pixels) of the digital image in Fig. [21 which 
represents the unique white hole. This is similar to considering the digital 
image as a cell comple^l^ 0] (see Fig. [21c). Unfortunately digital images 
are not 'clean', noise can create many unwanted holes which complicate the 
correct interpretation. One can ask for the edges we have to delete in order 



^Intuitively a cell complex is defined by a set of 0-cells (vertices) that bound a set of 
1-cells (edges), that bound a set of 2-cells (faces), etc. 
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to 'destroy' a hole. In the example in Fig. [21 it is not enough to delete only 
one edge. The deletion of the bold edges in Fig. |2Jc together with the faces 
that they bound produces the 'disappearing' of the hole. The set of bold 
edges in Fig. |2jc define a 1-cocycle, a topological invariant of the respective 
object. Equivalence classes of such cocycles are the elements of cohomology. 

To cope with complexity issues that arise when highly complex algorithms 
must be applied to huge amounts of data, graph pyramids are used. These 
hierarchical data structures offer possibilities to reduce the amount of data 
by local operations which can be applied in parallel and which have the 
enormous advantage to preserve the topological properties of the data. Hence 
the search for independent cocycles can be correctly done on a fraction of the 
data at the top of the pyramid. The simplified geometry of these cocycles 
can then be delineated top-down through the levels of the pyramid by again 
local processes up to the high accuracy of the base level. 

Maybe due to its more abstract nature, lacking a geometric meaning and 
due to its computation complexity, cohomology has not been yet widely ap- 
plied to pattern recognition and image processing. This paper is possibly 
the first attempt to use it in the context of digital images. For this pur- 
pose, in this paper we consider the best known environment which are 2D 
images whereas nD is the ultimate goal. Concepts related to cohomology 
can have associated interpretations in graph theory. Having these interpre- 
tations opens the door for applying classical efficient graph theory algorithms 
to compute and manipulate these features. Besides, for objects embedded in 
M^, homology- a wider used topological invariant, and cohomology groups 
are isomorphic. But the ring structure presented in cohomology character- 
izes the relations between 2-holes (cohomology generators of dimension 1), 
which homology does not. Indeed, dealing with homology and cohomology 
properties, representative cycles and cocycles, and their computation is quite 
different and doing this study in 2D gives important insights which should 
be relevant for extension to nD, n > 2. Initial results regarding this work 
have been presented in jl] . The current paper extends our earlier publication 
with detailed insights and proofs, and a refinement of the previous method 
that makes the obtained cocycles scanning and rotation invariant in the case 
of an identical discretization. 

The paper is organized as follows: Sections [2] and [3] recall graph pyramids 
and cohomology, and make initial connections. In Section |U preserving- 
topology properties in irregular graph pyramids are given. Section [5] presents 
the proposed method. Section [HI uses the properties of the proposed method 
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a) b) c) d) 

Figure 2: a) A 2D digital image /; b) its RAG; c) a cell complex associated to / (with 
bold edges, a representative cocycle); and d) the cell complex without the hole. 

to extend it and obtain scanning and rotation invariant cocycles. Section [7] 
concludes the paper. 

2. Irregular Graph Pyramids 

A graph (see for example [sl) is an ordered pair G = iV^E) comprising 
a set V of vertices and a set E of edges. Each edge e G -E is incident to 
two not necessarily distinct vertices v^w G V, written as e = {v,w). An 
edge e = {v,w) is said to be directed if the pair iy^w) is an ordered pair. 
An edge is said to be a self- loop ii v = w. The edge e = {v,w) is called a 
parallel edge iff 3e' E E,e' e s.t. e' = {v, w). A graph is called undirected 
if none if its edges is directed, and it is called planar if it can be drawn in a 
plane with no edges crossing (vertices are drawn as points and edges as lines 
connecting their incident vertices). Given a graph G = {V,E) removing an 
edge e E E will result in the graph G' = {V,E\ {e}), contracting the edge 
e = {v,w) implies removing it and identifying its incident vertices s.t any 
remaining edge previously incident to v or w is now incident to the unique 
vertex v = w. 

Given a decomposition of an object or image into regions a region ad- 
jacency graph (RAG) is an undirected graph that encodes the adjacency 
of regions in a partition. A vertex is associated to each region, vertices of 
neighboring regions are connected by an undirected edge. Classical RAGs 
do not contain any self-loops or parallel edges. An extended region adjacency 
graph (eRAG) is a RAG that contains the so-called pseudo edges, which are 
self-loops and parallel edges used to encode neighborhood relations to a cell 
completely enclosed by one or more other cells j^. The dual graph of an 
eRAG G is called a boundary graph (BG) and is denoted by G {G is said to 
be the primal graph of G). The edges of G represent the boundaries (bor- 
ders) of the regions encoded by G, and the vertices of G represent points 
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Figure 3: A digital image /, and boundary graphs Gg, Gio and G15 of the pyramid of /. 

where boundary segments meet. G and G are planar graphs. There is a 
one-to-one correspondence between the edges of G and the edges of (5, which 
induces a one-to-one correspondence between the vertices of G and the 2D 
cells (will be denoted by /acej§) of G. The dual of G is again G. The fol- 
lowing operations are equivalent: edge contraction in G with edge removal 
in (5, and edge removal in G with edge contraction in G. 

A (dual) irregular graph pyramid 0, a a a is a stack of successively 
reduced planar graphs P = {(Go, Gq), . . . , {Gn, Each level {Gk, Gk), 

< k < n, is obtained by first contracting edges in Gk-i (removal in Gk-i), 
if their end vertices have the same label (regions should be merged), and 
then removing edges in Gk~i (contraction in Gk~i) to simplify the structure. 
The contracted and removed edges are said to be contracted or removed 
(sometimes called removal edges) in (Gfc.i, Gfc_i). In each Gk-i and 
contracted edges form trees called contraction kernels. One vertex of each 
contraction kernel is called a surviving vertex and is considered to have 'sur- 
vived' to (GkiGk)- The vertices of a contraction kernel in level k — 1 form 
the reduction window of the respective surviving vertex v in level k. The re- 
ceptive field of V is the (connected) set of vertices from level that have been 
'merged' to v over levels 0, . . . ,k. The equivalent contraction kernel (ECK) 
of a vertex v is the tree obtained by replacing v and all is descendants with 
their corresponding contraction kernels. The vertices of the ECK of v form 
the receptive field of v. 

Algorithm!!] gives the main steps used to build a graph pyramid. In LineO 
the operations are performed on both graphs (adjacency and boundary) even 
if for simplicity only one of the graphs is mentioned in each step. Line |5lii 
removes self-loops bounding an empty face in the adjacency graph. Line|5liii 
simplifies the region boundaries in the boundary graph. 



^Not to be confused with the vertices of the dual of a RAG (sometimes also denoted 
by the term faces). 
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Algorithm 1 Build (dual) graph pyramid 

Input: image / /* pixels labeled 'object' or 'back- 

ground'*/ 

1: (Go,Go) = ((K),i?o),(^0,^o)) 

/*Vo associates a vertex to every pixel. Eq connects vertices correspond- 
ing to 4-connected pixels. Go and Go are dual.*/ 

2: k = 

3: repeat 

4: /*select edges to contract*/ 
T = 

for all V E Gk do 

i. select an edge {v,w) G Gk with v,w having the same label 

ii. T ^ T U {v, w) /*add edge*/ 
5: if T ^ then 

/* region merging, easier described in the adjacency graph:*/ 

i. {G', G') contract edges T of Gu (removal in Gk) 

/* simplification, easier described in the boundary graph:* / 

ii. (G", G") ^ contract pending trees in G' (removal in G') 

iii. {Gk+i,Gk+i) ^ contract one distinct edge incident to each ver- 

tex of degree 2 in G" (removal in G") 

6: k ^ k + 1 
7: until T = 

Output: Graph pyramid P = {(Go, Gq), . . . , {Gk^i, Gk-i)}- 



3. Homology, Cohomology and Integral Operators 



We refer to [10[ for an introduction to homology and cohomology. 

Intuitively, homology characterizes the holes of any dimension (i.e. con- 
nected components, 1-dimensional holes, etc.) of an n-dimensional object. It 
defines the concept of generators which, for example for 2D objects are sim- 
ilar to closed paths of edges surrounding holes. More general, /c- dimensional 



manifolds surrounding {k + l)-dimensional holes are generators [10[, and de- 
fine equivalence classes of (/c + l)-holes. Cohomology arises from the algebraic 
dualization of the construction of homology. It manipulates groups of homo- 
morphisms to define equivalence classes. Intuitively, cocycles (the invariants 
computed by cohomology), represent the sets of elements (e.g. edges) to be 
removed to destroy certain holes. See Fig. [2Jc for an example cocycle. 
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0- cells 

1- cells 

2- cells 
1-chain 



{vi,V2,V3,Vi} 
{61,62,63,64,65,66} 
{/l} 



homologous cycles 



1-boundary 



1-cycle 
1-cycle 



ei + 63 
dfi = 61 + 62 + 65 
a = 63 + 64 + 65 
6 = 61 + 62 + 63 + 64 
a and b; since a = b + dfi 



Figure 4: Example of cell, chain, boundary and cycle. 



A homeomorphism is a bijective continuous function between two spaces, 
that has a continuous inverse function. They are the mappings which preserve 
all the topological properties of a given space. Two spaces with a homeo- 
morphism between them are called homeomorphic, and from a topological 
viewpoint they are the same. 

Two continuous functions from one topological space to another are called 
homotopic if one can be 'continuously deformed' into the other. Two spaces X 
and Y are homotopy equivalent if there are maps f : X ^ Y and g : Y X 
such that gf is homotopic to idx and fg is homotopic to idy Observe 
that if two spaces are homeomorphic then they are homotopic. A homotopy 
invariant is a topological property which is invariant under homotopy. 

A p-dimensional cell (or p-cell, for short) is a topological space that is 
homeomorphic to the p-dimensional ball B^. A 0-cell is homeomorphic to 
a point, a 1-cell to an arc and a 2-cell to a disk. Roughly speaking, a cell- 
( or CW- ) complex is built by gluing together the basic building blocks called 
cells. 

Figures H] and [5] illustrate the following abstract concepts. 
3.1. Homology 

The notion of p- chain is defined as a formal sum of p-cells. The chains 
are considered over Z/2 coefficients i.e. a p-cell is either present in a p- 
chain (coefficient 1) or absent (coefficient 0) - any cell that appears twice 
vanishes. The set of p-chains form an abelian group called the p-chain group 
Cp. This group is generated by all the p-cells. The boundary operator is 
a set of homomorphisms {dp : Cp — > Cp_i}p>o connecting two consecutive 
dimensions. By linearity, the boundary of any p-chain is defined as the formal 
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0- CO chain 


{^1} 


1-cochain 


{61,64} 


1-coboundary 


S{vi} = {63,64} 


1-cocycle 


C = {61,62} 


1-cocycle 


d = {61,62,63,64} 


homologous cocycles 


c and d; since d = c + S{vi} 



Figure 5: Example of cochain, cocyclc and coboundary. 



sum of the boundaries of each p-cell that appears in the chain. The boundary 
of 0-cells (i.e. points) is always 0. A chain complex is the set of all the chain 

groups connected by the boundary operator: ■ ■ ■ Cp % Cp_i— )■ ■ ■ ■ — ^ 

Co^O. 

A p-chain a is called a p-cycle if dp{cr) = 0. If a = 9p+i(/i) for some 
{p + 1) -chain /i then a is called a p- boundary. Two p-cycles a and a' are 
homologous if there exists a p-boundary b such that a = a' + b. Denote the 
groups of p-cycles and p-boundaries by Zp and Bp respectively. For each p, 
dp-idp = 0. In other words, all p-boundaries are p-cycles {Bp C Zp). Define 
the p^^ homology group to be the quotient group Hp = Zp/Bp, for all p. 
Each element of Hp is a class obtained by adding each p-boundary to a given 
p-cycle a. Then a is a representative p-cycle of the homology class a + Bp. 

Since the chains are considered over Z/2 coefficients, the chain groups are 
vector spaces and the boundary operators are linear operators. The cycle and 
and boundary groups are just the kernel and image of such operator. The 
homology group is a quotien space. 

3.2. Cohomology 

Cohomology groups are constructed by turning chain groups into groups 
of homomorphisms and boundary operators into their dual homomorphisms. 
Define ap-cochain as a homomorphism e : Cp — ?■ Z/2. We can see ap-cochain 
as a binary mask of the set of p-cells: imagine you order all p-cells in the 
complex, (let's say we have n p-cells, and call this ordered set Sp). Then a 
p-cochain c is a binary mask of n values in {0, 1}", where n is the number 
of p-cells in the complex. When no confusion can arise, we will identify the 
p-cochain c with the set S of p-cells that are evaluated to 1 by c. 

The p-cochains form the set which is a group. The boundary operator 
defines a dual set of homomorphisms, the coboundary operator {6'^ : — )■ 
C^~^^}p>Q, such that 6^{c) = cdp^i for any p-cochain c. Since the coboundary 
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Figure 6: The cell complex K, K' and K"; and the homomorphisms (/)p,7rp,tp and 

Vpi "pi ''p' 

operator runs in a direction opposite to the boundary operator, it raises the 
dimension. Its kernel is the group of cocycles and its image is the group of 
cohoundaries. A p-cochain c is a p-cocycle if 5^c(= cdp^i) : Cp+i — )■ Z/2 
is the null homomorphism. A p-cochain ci is a p-coboundary if there exists 
a — l)-cochain e such that d = 6^~^e{= edp). Two p-cocycles c and c' 
are cohomologous if there exists a p-coboundary d such that c = d + d. 
The p*^ cohomology group is defined as the quotient of p-cocycle modulo p- 
coboundary groups, = Z^/B^, for all p. Each element of is a class 
obtained by adding each p-coboundary to a given p-cocycle c. Then c is a 
representative p-cocycle of the cohomology class c + B^, denoted by [c]. 

Definition 3.1. A set of p-cocycles {ci,...,c„} is a basis of representative 
p-cocycles if: 

• Any other p-cocycle c can he written as a linear combination of the 
p-cocycles of the set plus a p-coboundary, that is: 

n 

c = ^ XiCi + 6P-^e 

i=l 

where Aj = 0, 1 and e is a {p — l)-cochain. 

• None of the p-cocycles in the set can be written as a linear combination 
of the rest plus a p-coboundary (minimality) . 

Remark 3.1. //{ci, . . . , c„} is a basis of representative p-cocycles then {[ci], 
. . . , [cn]} is a basis of the p*'* cohomology group of the object. Each [q] is a 
cohomology generator. 
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Figure 7: The liomomorphisms ^LttIJ and (pp^T^i 



3.3. Integral Operator 

Starting from a chain complex, ■ ■ ■ Ci — t- Co 0, take a g-cell a 
and a (g + l)-chain a. An integral operator 11| is defined as the set of 
homomorphisms {(pp : Cp Cp+i}p>o such that 0q(cr) = a, (f)q{fi) = if /i is 
a g-cell different to a, and for all p ^ q and any p-cell 7 we have 0p(7) = 0. 
It is extended to all p-chains by linearity. 

An integral operator {0p : Cp — )■ Cp+i}p>o satisfies the chain-homotopy 
property iff 4>pdp+i(f)p = (f>p for each p. For 0p satisfying the chain-homotopy 
property, define Hp = idp + (pp-idp + 9p+i0p where {idp : Cp — )■ Cp}p>o is 
the identity. Define imiTp = {b E Cp \ 3a E Cp s.t. 7rp{a) = b}. Then, 

■ ■ ■ % irrnii % imiTn^ is a chain complex and {np : Cp — )■ irmip}p>Q is a 



chain equivalence [10|. Its chain-homotopy inverse is the inclusion map {tp : 
irrnip — )■ Cp}p>o. Integral operators satisfying the chain-homotopy property 
can be seen as a kind of inverse boundary operator: They raise the dimension 
and satisfy the nilpotent condition 0p+i0p = for all p. Although, in general, 
4>p-idp 7^ idp and Sp+if/'p 7^ idp, what happens is (j)pdp+i(f)p = (pp for all p 
(which would be equivalent to x ■ - ■ x = x for x G M \ {0}). Consider, for 
example, the cell complex K in Fig. E] on the left. The integral operator 
associated to the removal of the edge e is given by 0i(e) = A. Then, 
7r[(e) = {idi + <Pldi + 9200(e) = e + a + f + d + e = a + f + d, tiI{A) = 0, 
7r2(-B) = A {A + B is renamed as A in K') and vr^ is the identity over the 
other p-cells of i^, p = 0, 1, 2. The removal of edge e decreases the degree of 
vertex v allowing for further simplification. 

The following lemma guarantees the correctness of the down projection 
procedure for computing cocycles given in Section [51 Since graph pyramids 
offer possibilities to reduce the amount of data by local operations, then 
the search for independent cocycles can be done on the top of the pyramid. 
Hence, the following lemma guarantees that these cocycles can be correctly 
delineated top-down through the levels of the pyramid. 
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Lemma 3.2. Let {0p : Cp — > Cp+i}p>o be an integral operator satisfying the 

chain-homotopy property. The chain complexes ■ ■ ■ ^ Ci % Co % and 

■ ■ ■ % imiTi % imiTo have isomorphic homology and cohomology groups. 
If c : imTTp —7- Z/2 is a representative p-cocycle of a cohomology generator, 
then air : Cp ^ Z/2 is a representative p-cocycle of the same generator. 

Proof. An integral operator that satisfies the chain homotopy property, is a 
chain homotopy of the identity {idp : Cp — > Cp}p>o to {ipUp : Cp — > Cp}p>o. 
Therefore, {vTp : Cp — )■ im7rp}p>o is a chain equivalence and chain equivalences 
induce isomorphisms on homology and cohomology (see [loj)- D 
For example, consider the cell complex K' of Fig. El The 1-cochain a*, 
defined by the set {b, /} of edges of K', is a 1-cocycle which 'blocks' the white 
hole H (in the sense that all the cycles representing the hole must contain 
an odd number of edges of a*). Then (3 = a*Ti\ is defined by the set {6, /, e} 
of edges of K. a and /3 are both 1-cocycles representing the same white hole 
H. 

4. Preserving Topology in Irregular Graph Pyramids 

Considering binary images, an object is defined by a connected set (4- 
connectivity) of foreground pixels. A partition of the whole space (foreground 
and background) in cells is called a cell subdivision. The referred partition 
could be obtained from any of the planar graphs in every level of the pyramid. 

Fix a level (Gj, Gj), the cell complex associated to the foreground object, 
called boundary cell complex, denoted by i^*, is obtained from (Gj, Gi) by 
taking all faces of Gi corresponding to vertices of Gj, whose receptive fields 
contain (only) foreground pixels, and adding all edges and vertices needed 
to represent the faces. The p-chain group generated by the p-cells of is 
denoted by Cp{K'). 

The following lemma guarantees that the local operations applied to build 
a graph pyramid preserve the topological properties of the initial data. 

Lemma 4.1. The boundary cell complex is well-defined. All the boundary 
cell complexes of a given irregular dual graph pyramid are cell subdivisions of 
the same object. Therefore, all these cell complexes are homeomorphic. 

Proof. Our input is a binarized 2D digital image. An object is the 4- 
connected set of foreground pixels. Since we only remove an edge in Gk 



12 





4£J 



b) 



d) 




Figure 8: a) Boundary cell complex K4 obtained from the top level of the pyramid, 
(6*4,(54); b) in bold, spanning tree edges of 04; c) homology-generator level, K"; d) in 
bold, the self- loops representing the cocycle edges in K"; e) in bold, the cocycle edges in 
top level. 



when it is in the boundary of two different regions that have the same label 
(region merging), and contract an edge in Gk when it is incident to a vertex 
of degree 2 or it is a pendant edge (simplification), all the new p-cells created 
are homeomorphic to p-dimensional balls, p = 0, 1, 2. □ 

As a result of Lemma 14.11 topological invariants computed on different 
levels of the pyramid are equivalent. 

For the purpose of this paper, a new cell complex called homology-generator 
level, is added over the boundary cell complex obtained from the top (last) 
level of the pyramid (see Fig. [S]). This new cell complex is denoted by 
and it is a set of regions surrounded by a set of self-loops incident to a single 
vertex. To obtain this cell complex, on the top of the computed pyramid, 
we compute a spanning tree of the boundary graph of the top level of the 
pyramid, and contract all the edges that belong to it. Note that K'^ is no 
longer homeomorphic to any K\ but homotopic. 

Lemma 4.2. The two operations used to construct an irregular graph pyra- 
mid: edge removal and edge contraction, are integral operators satisfying the 
chain-homotopy property. 

Proof. Fix a level {Gi,Gi), suppose an edge e in Gi is removed. Since Gi 
is planar, then e is in the boundary of two 2-cells (or regions) A and B (see 
Fig. E]). The integral operator 0'' associated to this edge removal is given by 
(f>{{e) = B (see Fig. [7]). Now, suppose that an edge d of Gj, with a vertex v 
of degree 2 in its boundary, is contracted (see Fig. [6]). The integral operator 
(j)'^ associated to this edge contraction is given by ^qI"^) — ^ (see Fig. [7]). □ 
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Starting from a cell decomposition of an object, its homology studies 
incidence relations of its subdivision. Cohomology arises from the algebraic 
dualization of the construction of homology. Both homology and cohomology 
are homotopy invariants. 

Corollary 4.3. The boundary cell complex of any level of the pyramid and 
the homology- generator level have isomorphic homology and cohomology groups. 

As a consequence of Lemmas 13.21 and 14.21 we have: 

Lemma 4.4. Fix a level {Gi,Gi), suppose an edge e in Gi, which is in the 
boundary of a region B, is removed. Let 0i(e) = B be the integral operator 
associated to this removal. Let a* be a 1-cocycle defined by a set of edges S 
in Gi \ {e}. // an odd number of edges of a* is in B, then S U {e} defines a 
1-cocycle in ; otherwise, it is S which defines a 1-cocycle in . 

In terms of embedded graphs, an integral operator maps a vertex/point 
to exactly one of its incident edges and an edge to exactly one of its incident 
faces. In every level of a graph pyramid, the contraction kernels make up 
a spanning forest. A forest composed of k connected components, spanning 
a graph with n vertices, has k root vertices, n — k other vertices, and also 
n — k edges. These edges can be oriented toward the respective root such that 
each edge has a unique starting vertex. Then, integral operators mapping 
the starting vertices to the corresponding edge of the spanning forest can be 
defined as follows: (poivi) = cj, where Cj is the edge incident to Vi, oriented 
away from it. 

The following lemma guarantees that all integral operators that create 
homeomorphisms are in fact a combination of the two operations used to 
construct irregular graph pyramids. 

Lemma 4.5. All integral operators that create homeomorphisms can be rep- 
resented in a dual graph pyramid. This is equivalent to: given an input 
image {Go,Go) and its associated cell complex Z = {Go,Gi,G2}, a cell com- 
plex Z' = {Cq,G[,G2} with Z,Z' homeomorphic, and Z a refinement of Z' 
i.e. Cq C Co, G[ C Ci, and G2 C G2, then there exists a pyramid P s.t. Z' 
is the cell complex associated to some level {Gk,Gk), k > 0, of P. 
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Level 



Level 1 Level 2 Level 3 Level 4 










Figure 9: Top row, from left to right: boundary graphs for all levels of the pyramid. 
Vertices surviving to the next level are drawn with a square. Middle row. in bold: con- 
tracted edges in the respective levels. Bottom row: removed edges in the boundary graph, 
equivalent to contracted edges in the adjacency graph. 



5. Representative Cocycles in Irregular Graph Pyramids 

A method for efficiently computing representative cycles of homology gen- 



erators using an irregular graph pyramid is given in [12| . In [13[ a novel algo- 
rithm for correctly visualizing graph pyramids, including multiple edges and 
self-loops is given. This algorithm preserves the geometry and the topology 
of the original image and has been used to produce the images throughout 
the paper (see Fig. ^ 

In this paper, representative cocycles are computed and drawn in the 
boundary graph of any level of a given irregular graph pyramid. They are 
computed in the top level and down projected using the described process. 



In the homology-generator level (see Fig. Elc), each self- loop a that sur- 
rounds a region of the background (hole of a region R of the foreground) is 
a representative 1-cycle of a homology generator. 
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Let be the homology-generator level. Without loss of generality, we 
can suppose that K" is connected. If not, repeat the following reasoning for 
each connected component (region) of K" . Let {ai, . . . , an} be the set of the 
self-loops surrounding a face of the background. Therefore, there are n white 
holes: Oi, . . . (see Fig. [Hid). Fix z, z = 1, . . . , ra, aj is a representative 1- 
cycle of the homology generator associated to the white hole Oj. Let /3 be a 
self-loop surrounding the face / of the foreground (recall that we suppose that 
K" is connected) such that ai is in the boundary of / in K" . Form the sets 
{«!, /?}, . . . , Let denote the boundary cell complex associated to 

the foreground in Gq. Let {(pp : Cp{K^) Cp+i{K^)}p>o be the composition 
of all integral operators associated with all removals and contractions of 
edges of the foreground of the boundary graphs of a given irregular graph 
pyramid. Let {np = idp + (pp-idp + 9p+i</)p : Cp(A'°) — )■ Cp{K")}p>o where 
{ip : Cp{K") — > Cp{K^)}p>Q is the inclusion map. 

Proposition 5.1. The 1-cochain a* defined by the set in is a 1- 

cocycle. Moreover, the set {aj, . . . , a*} is a basis of representative 1-cocycles. 

Proof. The set of edges of K" is the set of the self-loops {ai, . . . , a„} sur- 
rounding a region of the background together with the self- loop f3, renamed 
by On+i, surrounding the face /. 

The 1-cochain a* is a cocycle in K" iff Si{a*) = 0. Since we work with 
objects embedded in then ai can only be in the boundary of two faces. In 
this case, one face belongs to the background and the other face is / in K" . 
Then 5i («*)(/) = a*id,if)) = «*(a, + + ■ ■ ■) = a*(ai) +<(/?)+ <(■ " ■) = 
1 + 1 + = 0. 

Let us prove minimality. Suppose, for example, that al = a*^ + ■ ■ ■ + a*^ 
where 1 < ji < ... < Js n, s > 1. Then a J is defined by the set 
{oj^, . . . , ftj^} if s even, and {aj^, . . . ,aj^, f3} if s odd, which is a contradiction. 

□ 

We will say that a* is a representative 1-cocycle of the cohomology gen- 
erator associated to the white hole Oj. 

Algorithm [2] gives the proposed method to downproject a cocycle a* from 
level k to level k — 1. Informally, in the homology-generator level, there is 
only one face representing the object. Based on the geometric interpretation 
of cocycles (Section [3]), if we remove the edges and the face in-between, we 
destroy the hole. Then, there is no need to add any other edge to the cocycle 
to remove the hole. However, when going down in the pyramid, this face 
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Algorithm 2 Down project cocycle 

For each connected component (region) of K" . Let A^, k > 0, denote the 
set of edges that define a cocycle in Gk (the boundary graph in level k). 
The down projection of Ak to Gk-i is the set of edges A^-i C Gk-i that 
corresponds to A^ i.e. represents the same cocycle. A^-i is computed as 
Ak_i = Al_-^ U where Al_^ denotes the set of surviving edges in Gk-i 

that correspond to Ak, and v4^„^ is a subset of removed edges in Gk-i- The 
following steps show how to obtain Al._^: 

1. Consider the contraction kernels of Gk-i (RAG) whose vertices are 
labeled with i (the region for which cocycles are computed). The edges 
of each contraction kernel are oriented toward the respective root - each 
edge has a unique starting vertex. 

2. For each contraction kernel T, from the leaves of T to the root, let e be 
an edge of T, v its starting point, and Ey the edges in the boundary of 
the face associated to v: label v with the sum of the number of edges 
that are in both Al._^ and in Ey, plus the labels of the children nodes 
of V. 

3. A removal edge of Gk-i is in if the starting point of the corre- 
sponding edge of Gk-i is labeled with an odd number. 



is partitioned. A connection among all the new regions is determined by 
the contraction kernels of the RAG. When the first partitioning occurs, the 
contraction kernel will contain one or two nodes corresponding to faces with 
one surviving cocycle edge in its boundary, and the rest of the nodes will have 
none. What Algorithm [2] does is to find the unique path in the contraction 
kernel joining these two nodes, and take the set of boundary edges between 
consecutive faces as part of the new cocycle (see Prop. \6.S\ Ig.gj and their 
proofs). Lower levels will update the connections in subsections of the cocycle 
path. Every subsection will correspond to a partitioned region between two 
consecutive cocycle edges. 

Consider the example in Fig. El In the homology-generator level we have 
A5 = {a, (3} the representative 1-cocycle of a cohomology generator (self- 
loops in Fig. (Hid). For down projection in level 4, A4 = A^U A^. We have 
that Al is the surviving edges in bold of top level in Fig. |8le). In this case, 

= because there is no merging of foreground regions from the boundary 
cell complex obtained from the top level to the homology-generator level. 
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Figure 10: Top row: removed edges in boundary graph. Bottom row, from right to left: 
the down-projected cocycles in bold. Filled circles on faces, represent surviving vertices 
from the adjacency graph in foreground regions. 



In the example in Fig. [TOl the cocycle a* in level 4 is the set of the two 
edges in bold (see Fig. [101 bottom row, column d). The down projection 
from level 4 to 3 are the surviving edges of the cocycle in level 4. This is 
because there was no contraction in the foreground region. The contractions 
of the adjacency graph can be seen in the top row of the figure. 

In level 2 (Fig. [TOl bottom row, column c), the first contraction of fore- 
ground in the adjacency graph with a single edge appears. In this case, 
the leaf node represents a face with an even number (2) of surviving cocy- 
cle edges in its boundary, which leads to not adding any other edge to the 
down-projected cocycle. Only in the base level (column a), one contraction 
kernel has a leaf node with an odd number (1) of surviving cocycle edges in 
the boundary. In this case, the corresponding edge in the boundary graph, 
for the edge connecting the respective node with its father, is added to the 
cocycle. 

Any edge Gq that has survived to a higher level and was selected as 
part of the cocycle in Gk-, will belong to the down-projected cocycle in Gq. 
In particular, the edges in Go that have survived to be the edges a, /3 of the 
cocycle in the homology-generator level, are going to be the entry and exit 
point of the cocycle path through the foreground region. 

In Fig. [m the space between the outside boundary of the object and the 
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Figure 11: a) Boundary cell complex obtained from the top level of the pyramid, (G4, G4); 
b) homology-generator level with cocycles edges in bold; c) ECK of the foreground in 
RAG, drawn over Go with cocycle edges in bold; d) to g) shows removed edges in bold for 
levels from 3 to 0. Cocycle edges arc marked with two small parallel lines. 



hole is bigger, allowing for more possibilities for the paths of the cocycle. The 
cocycle path in the base is going to converge to the unique path connecting 
the surviving edges a, (3 through the ECK (see Fig. [TTlc). 

In Fig. [Hid the cocycle is made of the surviving cocycle edges from the 
homology-generator level in Fig. [TTlb . The first partition of the foreground 
is connected by the contraction kernel in Fig. [11] e . Here, one of the regions 
in the partition contains in its boundary the two surviving cocycle edges, so 
there is no path to find and no new cocycle edge to add. In Fig. [TTlf there 
is only one edge to add to the cocycle to connect the path of edges, which is 
identified with the leaf node with label 1. 

Notice that from level Fig. [TTJf to the one in Fig. [Hlg also the selected 
'surviving' edges play a role, as edges to be removed at higher levels had to 
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have been 'surviving' ones at levels below. In this case, the edge we add to 
the cocycle in level 1 (Fig. [Hlf), was a surviving edge in level below that was 
connecting two contraction kernels. Therefore, the selection of the surviving 
edges also determine the delineation of the down-projected cocycle together 
with the contracted ones. The ECK contains the decision of which edges 
were contracted but also which ones were simplified, determining the unique 
path. 

Proposition 5.2. 1. The down projection of a* is a set of edges 'blocking' 
the creation of the hole Oi, i.e., given a cycle g homologous to the down 
projection of the cycle ai, the down projection of a* contains an odd 
number of edges of g. 

2. The down projection of a* is always a cocycle. Moreover, the down 
projection of {a^, . . . , a*} is a basis of representative 1-cocycles. 



1. The down projection of a*, which is a*7ri, contains an odd number of 
edges of g iff «*((?) = 1. First, if g is homologous to the down projection 
of Qfj, which is ti(aj), then there exists a 2-chain b in Kq such that 
g = Li{ai) + ^2(6). Second, a*7ii{g) = a*7ii{Li{ai) + di{b)) = 1, since 
a*Li{ai) = = 1, and a*'7iid2{b) = because a* is a cocycle and 

ijjd2 = i9ivr2 (since {tt^ : Cp{K^) — j- Cp{K*H)}p>Q is a chain equivalence 
|lO|). So g must contain an odd number of edges of the set that defines 




2. Proof of correctness of the down projection algorithm: it is a conse- 



Example down projections are shown in Fig. |9l [TOl and [TT] 
5.1. Complexity 

Let n be the height of the pyramid (number of levels) and Vq, Cq the 
number of vertices, respectively edges in the base level, with n ^ logfo 
(logarithmic height). An upper bound for the computation complexity is: 
0{vQn) to build the pyramid; for each foreground component, 0{h) in the 
number of holes h to choose the representative cocycles in the top level; 
O(eon) to down project each cocycle. The overall computation complexity 



Proof. 




a*. 



quence of Lemma 14.41 



□ 
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a) large hole - outer boundary b) small hole - outer boundary 



Figure 12: Example showing (top) normal and (bottom) rotation invariant eocycles. The 
cocycles are down-projected starting with the surviving edge of the outer boundary and 
the surviving edge of the a) large and b) small holes. 

is then below 0{von + c{heon)), where c is the number of cocycles that are 
computed and down projected. 

Actually not all edges are part of cocycles and not all levels have cq 
number of edges. When building Ak-i one can go in linear time over the 
edges of and consider only the contraction kernels in Gk-i for which the 
surviving vertices in correspond to one of the two faces to which an edge 
of Ak is incident to. Then, computing Ak_i actually takes \Ak\ + l^jt-il 
number of steps, where T^„i G Gk^i are the contraction kernels mentioned 
before. Thus in practice the complexity of down projecting a cocycle is below 
O(eon). 

6. A First Step Towards Stable Cocycles 

If topology is considered in the context of recognition, or a joint extraction 
of both topological and geometrical features is required, then the location 
and shape of the extracted topological invariants becomes relevant and is an 
important way to ensure stability/repeatability. A relevant example is the 
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work in [IJ] where handle and tunnel loops are made geometry aware by 
placing them on "geometrically relevant" positions. 

In this section we make a first stept towards obtaining stable cocycles 
and consider invariance with respect to scanning and rotation of the object. 
Because we start with a pixel grid and the 4 neighborhood, only rotations with 
multiple of 90 degrees produce identical discretizations allowing for identical 
cocycles. 

The following properties are required for Prop. 16.41 which gives the parts 
of the pyramid that the computed cocycles depend on. 

Proposition 6.1. The down-projected cocycles contain only removal edges 
in the boundary graphs, corresponding to edges in the ECK of the top vertex 
representing the object in Gn, and the two edges a, /3 that have survived to 
Gn and where selected as the cocycle a* = {a,f3}. 

Proof. Algorithm |2] starts with two edges {a, (3} in Gn, and for each level 
/c = n — l,...,Oit adds only removal edges from Gk i.e. edges that where 
contracted in Gk to merge neighboring regions belonging to the object. □ 

Proposition 6.2. The result of down projection (Algorithmic does not de- 
pend on the selected surviving vertices in Gk- 

Proof. Consider the function q : T C G^^i N, q{T) = J^vgt l^fc-i ^ ^^1 
(see Alg|2]for the used notation). Every cocycle Aj. has an even number of 
edges from the boundary of any face in Gk (Section [3]). Then q{T) is also even 
and the number of vertices v for which \Al_-^ fl E^l is odd, is even. For any 
edge e G T = iV,E) consider the two connected components (trees) Ti,T2 
of the subgraph (V, E \ {e}) (e is a cut edge of T because T is a tree). The 
removal edge of Gk^i corresponding to the edge e is added to the cocycle 
if q{Ti) and q{T2) are even, which is independent of the originally chosen 
surviving vertex, the root of T. □ 

Proposition 6.3. A down-projected cocycle does not depend on the order in 
which edges are removed in the boundary graphs (region merging). 

Proof. Consider the proof of Prop. 16.21 What Algorithm [2] does is to select 
additional edges from the contraction kernel T, to connect the vertices v 
with odd \Al_-^ fl E^\. In a tree, there is a unique path connecting any two 
vertices. Denote by K the ECK of the vertex in G„ corresponding to the 
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a) 



b) 



Figure 13: a) original image; b) in bold, the paths in the RAG Go associated to the down- 
projected cocycles related to the holes representing the top-left and top-right windows. 



object. Algorithm|2]returns a and /3, plus the set of edges of Gq corresponding 
to the path in K that connects the two vertices whose corresponding faces in 
Go have a and /3 in their boundary. The ECK of a vertex does not depend 
on the order of the intermediate steps 15||. □ 
The following property is an immediate result of Properties 16.11 16. 2[ 
and l6.3[ 



Proposition 6.4. The cocycles computed by Algorithmic depend only on the 
cocycle {a, f3} chosen in the top level, {Gn, Gn), and on the ECK of the vertex 
in Gn corresponding to the face describing the object in Gn- 

The following property results from Prop. 16.41 and motivates the modifi- 
cation proposed in the rest of this section. 

Proposition 6.5. If the ECK of the vertex representing the object in the ad- 
jacency graph of the homology- generator level, and the edges that survive to 
he in the boundary of the corresponding face in the boundary graph, are scan- 
ning and rotation invariant, we will obtain scanning and rotation invariant 
cocycles. 



In the following we will consider the necessary additions to the pyramid 
building process, to ensure that computed cocycles do not depend on the 
scanning and rotation of the object. We follow Prop. 16.51 and consider the 
ECK and the boundary edges of the top level. 
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b) 



Figure 14: a) original image; b) in bold, the paths in the RAG Go associated to the 
down-projected cocycles. 

6.1. Invariant ECK and surviving edges 

As the edges to be removed in Gk are 'locally' chosen, in a binary image 
like the ones used to represent our objects, there is no local structure and 
a random or scanning/ orientation dependent direction is taken (Line H] of 
Algorithm [1]) . 

To 'add structure' and create an ordering for selecting edges to be removed 
and contracted: 

1. Compute a spanning tree of the subgraph O C Go corresponding to 
our object in the base level. Mark the edges of Gq corresponding to the 
edges of the spanning tree to he the removed ones. 

2. Create a strict ordering between any two edges. This ordering is used 
to select surviving edges during simplification, and thus controls the 
choice of edges in the homology-generator level. 

The tree. Given the graph O corresponding to our object, and a vertex s G O, 
we define d{v) : O — j- N to be the number of edges of the shortest path 
connecting v and s in O i.e. the geodesic distance between the two pixels 
corresponding to v and s, using the 4 neighborhood. Note that a vertex 
s can be obtained in a rotation invariant manner for example by using an 
automatic shape orientation method 16| and then selecting the top, left-most 
vertex. 
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Stable ECK (in Gq). Every vertex v G Gq, f 7^ s, labels the edge (t>, v') G Gq 
with d{v') = d{v) — 1 as 'to contract'. If v has more than one neighbor v' 
with d{v') = d{v) — 1 the neighbor that minimizes the angle SVV , and in 
case of angle equality, the one that has a clockwise orientation of SVV is 
chosen. 5, V , and V are the points in 1? corresponding to the centers of 
the pixels represented by s, f , and v' . The edges in Go corresponding to the 
edges of Go labeled as 'to contract', are marked as 'to remove'. 

Region boundary simplification (in Gk)- In Line Oiii) of Algorithm [H from 
any chain of edges bounded by at least one vertex of degree 2, one edge 
will survive and all others will be contracted. To choose to surviving edge, 
assign to each edge e G Go bounding a cell of the object, the value /(e) = 
min{d{vi) , d{v2)} where vi, V2 are the vertices of Gq corresponding to the 
two faces of Go to which e is incident to. Faces not part of the object are 
ignored. When choosing the edge to survive i.e. not contract, the edges are 
sorted using the following (transitive) relation between any two edges e and 
e': 



. /(e) vs. fie'); 

• if /(e) = /(e') then use the orientation of c(e) S c{e') vs. the orientation 
of c(e') S c{e), where c(e), c(e') G are the centers of the edges e and 
e', and S* G is the center of the pixel used to define the rotation 
invariant tree; 

• if c(e) S c(e') are collinear, the Euclidean distance between c(e) and S 
vs. the Euclidean distance between c(e') and S. 

Homology- generator level. When building the homology-generator level, all 
edges of G„_i bounding the face of the object are sorted based on the criteria 
above. Edges are selected in inverse order and used to create the spanning 
tree to be contracted. Edges not bounding the face corresponding to the 
object are added in random order. 

Fig. [12] shows an example object and its computed cocycles with and 
without the rotation invariant pyramid. Fig. [T3] and [T3] show the paths in 



the RAG associated to the down-projected cocycles of the test images in [12 
Fig. US] is another example showing the path in the RAG Go associated to the 
down-projected cocycle related to the hole associated to the ball. Finally, 
Fig US] shows the down- projected cocycles computed on a image from the 
2010 World Cup final in South Africa. 
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Figure 15: In bold, the path in the RAG Go associated to the down-projected cocycle. 



6.2. Discussion 

Besides stability with respect to geometric transformations, one could 
consider additional criteria like the minimality of the obtained cocycles w.r.t. 
additional measures, like for example the number of edges of each cocycle, 
the sum of the number of edges of a basis of representative cocycles, the 
length of the path in passing through the support squares of the pixels 
having at least one cocycle in their boundarjjf], etc. 

It has been shown for homology generators 17| that in general the prob- 
lem of computing minimal representative cycles is NP-hard. Nevertheless, in 
certain cases, like computing (n — l)-cycles for n dimensional objects, finding 
minimal cycles is not NP-hard. For cocycles such a study does not exist 
yet, but considering the relation between homology and cohomology, similar 
results can be expected. In the case of 2D objects, the problem of finding co- 
cycles with minimal number of edges can be related to the problem of finding 
in the RAG shortest paths that connect vertices adjacent to different holes 
- this problem can be solved in nlogn if using Dijkstra's algorithm 



Real life objects are typically obtained by using scanning devices of dif- 
ferent type: 3D scanners, video cameras, CT, MRI, etc. One common issue 



■^Would give Euchdean rotation robust cocycles. 
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in all these cases is the presence of noise and being robust w.r.t. the possible 
deformations of a real object. A possible solution could be to define a ro- 



bust basis of cocycles based on a function like the eccentricity transform [19 
which is known to be robust w.r.t. noise and deformations. The eccentricity 
transform associates to each point of a shape the geodesic distance to the 
point furthest away. For a given starting point the geodesic distance function 
defines behind holes a set of points called the cut locus which can be reached 
in the same distance on multiple paths (going on both sides of the hole). In 
many cases, cutting a shape along these sets can produce a shape with less 
holes, which gives the same geodesic distance function for the same starting 
point. The eccentricity transform can be interpreted as the maximum over 
multiple geodesic distance propagations initiated at each point of the shape. 
As cocycles can be seen as 'cuts into the object' that 'kill' a hole, a cut 
which has the minimum effect on the eccentricity transform, could provide 
an avenue for selecting robust cocycles. 

Finding associations between concepts in cohomology and graph theory 
will open the door for applying existing efficient algorithms (e.g. shortest 
path). The following lemma can be seen as a first step in this direction. 

Lemma 6.6. Any set of foreground edges in the boundary graph Go, asso- 
ciated to a path in the RAG Go, connecting a hole Oi of the object with the 
(outside) background face, is a representative 1-cocycle cohomologous to the 
down projection of the 1-cocycle a*. It blocks any generator that would sur- 
round the hole. 

In other words, consider the down projection JlA] of Ui and (3 in Gq: the 1- 
cycles Li{ai) = a and ii(/3) = b, respectively. Take any edge Ca E a and Cb G b. 
Let fa, fb be faces of , the boundary cell complex associated to the fore- 
ground in Go having Ca, respectively Cb, in their boundary. Let vq, fi, . . . , Vn 
be a simple path of vertices in Go s.t. all vertices are labeled as foreground, 
vo is the vertex associated to fa, and f„ to fb. Consider the set of edges 
c = {co, . . . , Cn+i} of Go, where cq = Ca, e„+i = Cb, and ei, i = 1 . . .n, is the 
common edge of the regions in Gq associated with the vertices and Vi. c 
defines a 1-cocycle cohomologous to the down projection of the 1-cocycle a*. 

Proof, c is a 1-cocycle iff 082 is the null homomorphism. First, cd2{fe) = 
c{e£ + e^+i) = 1 + 1 = 0. Second, if / is a 2-cell of K^, f fg, i = 0, . . . ,n, 
then, cd{f) = 0. To prove that the cocycles c and a*7ri (the down projection 
of a* to the base level of the pyramid) are cohomologous, is equivalent to 
prove that cli = a*. We have that cti(a) = c(eb) = 1 and cti(/3) = c(ea) = 
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Figure 16: On the top-left, original image (a image from the 2010 World Cup final in 
South Africa); in bold, down-projected cocyclcs. 

1. Finally, cli over the remaining self- loops of the boundary graph of the 
homology-generator level is null. Therefore, CLi = a*. □ 

7. Conclusion 

This paper considers cohomology in the context of graph pyramids. Rep- 
resentative cocycles are computed at the reduced top level and down pro- 
jected to the base level corresponding to the original image. Connections 
between cohomology and graph theory are proposed, considering the appli- 
cation of cohomology in the context of classification and recognition. The 
current paper extends the previous work with detailed insights and proofs, 
and a refinement of the previous method that makes the obtained cocy- 
cles scanning and rotation invariant. Extension to higher dimensions, where 
cohomology has a richer algebraic structure than homology, and complete 
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cohomology - graph theory associations are proposed for future work. For 
this last task, we could consider nD generalized map pyramids 20 . 
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